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Abstract 

We present the measurements of the higher moments of the net-proton multiplicity distributions 
for Au+Au collisions from the first phase of the beam energy scan program at RHIC. The mea- 
surements are carried out at midrapidity {\y\ < 0.5) and within the transverse momentum range 
0.4 < pt < 0.8 GeV/c by STAR. The moment products , kct^ and Scr, which are related to the 
ratios of various order baryon number susceptibilities and sensitive to the correlation length , are 
applied to search for the QCD critical point. We report the centrality and energy dependence of 
the moment products (kct^ and Scr), results are compared with Poisson expectations. The ex- 
periment data are also compared with transport model calculation to understand the non-critical 
physics effects. 



1. Introduction 

In the year 2010, the Relativistic Heavy Ion Collider (RHIC) launched the Beam Energy 
Scan (BES) program |T] and tuned the Au+Au collision energy from y =200 GeV down to 
7.7 GeV. The main goal of the BES program is to look for the signature of the first order phase 
transition and Quantum Chromodynamics (QCD) Critical Point (CP) 12. S B S • By varying 
the colliding energy, we can access a broad region of the QCD phase diagram with a baryon 
chemical potential (/v^) range from about 20 MeV to 450 MeV [JJ. As the diverges of the 
correlation length (^) near the CP, non-monotonic variations of observables with collision energy 
are believed to be the best signature of CP. Higher moments (skewness {S ) and kurtosis (/c) etc.) 
of conserved quantities, such as net-baryon, net-charge and net-strangeness, distributions are 
more sensitive to the correlation length la, Q. Si than the variance (cr^) and could be detectable 
if the signature survive in the evolution of the system. On the other hand, the moment products, 
Ko^ and S cr, are also related to the ratios of various order susceptibilities as in Lattice QCD |@, 
Tol n], HRG model lfl2ll and PQM model ifisll . Ratios of baryon number susceptibilities used 



to compare with the experimental data as kct^ - x^b'IXb' Scr - X'b'I^^'b cancel out the 
volume effect. Theoretical calculations demonstrate that the experimentally measurable net- 
proton (proton number minus anti-proton number) number fluctuations can effectively reflect the 
fluctuations of the net-baryon number fTT). In this paper, we will present measurements for 
higher moments of event-by-event net-proton multiplicity distributions from the first phase of 
flie BES program at RHIC. 
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2. Results and Discussion 



The data of Au+Au collisions at ^/i;;;^!.!, 1 1.5, 19.6, 27, 39, 62.4 and 200 GeV from the 
first phase of the BES program are used in the analysis. We have also studied p + p collisions 
at -y/i^=62.4 and 200 GeV. The protons and anti-protons are identified with the ionization en- 
ergy loss (dE/dx) measured by Time Projection Chamber (TPC) at midrapiidty {\y\ < 0.5) and 
within the transverse momentum range 0.4 < pj < 0.8 GeV/c. To suppress autocorrelation 
effects between measured net-proton fluctuations and centrality in our results, centralities are de- 
termined by the uncorrected charge particle multiplicity (dNch/drj) by excluding the protons and 
anti-protons within pseudorapdity \ri\ < 0.5, a new technique add to the moment analysis after 
QM2011 conference. If protons and anti -protons are independent Poissonian distributions, the 
net-proton multiplicity will follow the skellam distribution, which is expressed as: 

^(A^) = (.^f'^lN(2 ^[M;MJ) exp[-(Mp + Mj)], 

where is a modified Bessel function, Mp and Mp are the measured mean values of protons 
and anti-protons. The various order cumulants (C„) are closely connected with the moments, 
e.g., Ci N MXi =< (AA^)^ >= cP-,Ci =< (AA^)^ >= ScP,Ca i^Nf > -3 < 
(AN)^ >^= Kcr'^. If the net-proton follows the skellam distribution, then we have, Scr = C3/C2 - 
(Mp - Mp)/{Mp + M-p) and kcP- - C4/C2 = 1, which provide the Poisson expectations for the 
moment products. 




Figure 1 : (Color Online) Centrality dependence of the cumulants of net- 
proton distributions for Au+Au collisions at yfT^=l.l, 11.5, 19.6, 27, 
39, 62.4 and 200 GeV. The lines are linear fit to these cumulants. Errors 
bars are statistical and caps are systematic errors. 



Figure 2: (Color Online) Centrality dependence of Scr/Poisson and KfP' 
of net-proton distributions for Au+Au collisions at y.v^^=7.7, 11.5, 
19.6, 27, 39, 62.4 and 200 GeV. The error bars are statistical and caps 
are systematic errors. 



The various order cumulants (Ci - C4) shown in Fig. [T]are obtained from the net-proton 
distributions and correct for the finite centrality bin width effect [4J. The statistical error estima- 
tion in Figlllare based on the Delta theorem method lITsll . which is an updated method compare 
to the analytical method used in QM2011 results The systematical errors are estimated by 
varying the track quality cut condition and particle identification. C\ — C4 values in FiglT]show 
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linear increase with Np-m in each centrality. The lines shown in the FiglUare linear fit to those 
cumulants, which is motivated by the Central Limit Theorem (CLT). In Fig. [T] it is observed 
that the odd (Ci and C3) and even (C2 and C4) order cumulants group together. The difference 
between the two groups are smaller at lower energies. The ratios of the cumulants, which are 
connected to the moment products as Scr - C3/C2 and Kcr^ = C4/C2, are shown in Fig. [S] It 
is observed that the kct^ and the S cr values normalized to Poisson expectations are below unity 
for y[s^ above 11. 5 GeV and above unity for 7.7 GeV in Au+Au collisions. The data presented 
here also provide information to extract freeze-out conditions in heavy-ion collisions using QCD 
based approaches Olol . 
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Figure 3: (Color online) Energy dependence of k(t and Scr for net-proton distributions for four collision centralities 
(0-5%, 5-10%, 30-40% and 70-80%) measured at STAR. The results are compared to UrQMD model calculations and 
p -V p collisions at yfT^=62A and 200 GeV. The lines in top panel are the Poisson expectations and in the bottom panel 
shows the 5(T normahzed to the corresponding Poisson expectations. 

Fig[3] shows the energy dependence of kcP- and ^cr of net-proton distributions for four cen- 
traUties (0-5%, 5-10%, 30-40% and 70-80%) in Au-hAu collisions. Bottom panel of Figl3]shows 
Scr values normalized to the corresponding Poisson expectations. The kct^ and normalized Scr 
values are similar for all of the centralities in Au-i-Au collisions at -^5^=39, 62.4 and 200 GeV 
and differences are observed between the 0-5% central Au-i-Au collisions and the peripheral col- 
lisions below -\fs^=39 GeV. The peripheral collision results are closer to unity and for -\fs^-7.7 
GeV, those are above unity. The UrQMD model iITtIi results are also shown in the Fig. [3]for 0-5% 
centrality to understand the non-CP effects, such as baryon number conservation and hadronic 
scattering. The UrQMD calculations show a monotonic decrease with decreasing beam energy. 
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3. Summary 

We present the beam energy ( ^JT^=^ .7 -200 GeV) and centrality dependence for the higher 
moments of net-proton distributions in Au+Au colUsions from the first phase of the BES pro- 
gram at RHIC. We have discussed the application of those observables to search for the QCD 
critical point. It is observed that the moment products, kct^ and S cr, are below the Poisson ex- 
pectation for central Au-hAu colhsions. The peripheral collision results are closer to unity and 
for GeV, those are above unity. The experimental data are also compared with the 

UrQMD model calculation to understand the non-critical physics effects. The UrQMD calcu- 
lations show a monotonic decrease with decreasing beam energy. Due to the small correlation 
length and dynamical evolution of the system in the heavy ion collisions, the signature of the CP 
could be small. This needs comparison between our data and QCD calculations with dynamical 
evolution of the system. We also need more statistics to get precise measurements below 19.6 
GeV and additional data at ■\/s^=l5 GeV. These are planned for the second phase of the BES 
program at RHIC. 
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